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SUMMARY 

Methods  have  been  developed  for  predicting  the  blast  loads,  quasi-static  pressure  rise, 
and  duration  of  the  quasi-static  pressure  rise  within  a suppressive  structure.  In  addition,  the 
side-on  overpressures  outside  a vented  suppressive  structure  are  estimated. 

Before  any  prediction  procedure  could  be  developed,  the  effective  vented  area  ratio 
(Oj,,  = - IveniM ) for  a multiwalled  structure  with  various  size  vents  in  each  wall  had  to  be 
developed-  This  relationship  was  assumed,  and  then  employed  in  a model  analysis  to  develop 
pi  terms  for  predicting  loads  and  durations  inside  the  structure  as  well  as  overpressures 
outside  the  structure.  Finally,  experimental  test  data  taken  from  the  literature  were  used 
to  develop  functional  relationships.  Wliereas  in  the  past,  investigators  have  assumed  that  the 
influence  of  a vented  suppressive  structure  was  a reduction  in  effective  cliarge  weight,  this 
solution  shows  that  a more  accurate  concept  is  the  creation  of  an  effective  standoff  distance 
less  than  the  frec-field  standoff  distance  at  which  blast  pressures  are  the  same  for  a given 
size  energy  release.  This  effective  standoff  distance  is  a function  of  the  effective  vented  area 
ratio  Oeii.  the  free-field  standoff  distance  fora  given  overpressure,  and  the  width  of  a 
suppressive  cubical  structure.  Data  from  a variety  of  test  structures  show  that  the  procedure 
predicts  outside  pressures  to  within  one  standard  deviation  of  18.6^ 

For  the  vast  majority  of  suppressive  structures,  the  quasi-static  pressure  rise  within  the 
structure  is  independent  of  the  vented  areas.  Test  data  indicate  that  provided  (Oeffvtj^'VF  is 
less  than  0.0775.  the  maximum  internal  quasi-static  pressure  is  a function  only  of  the  charge 
weight  divided  by  the  volume  because  the  maximum  pressure  is  reached  before  insignificant 
venting  occurs. 

This  report  is  a reprint  of  a paper  presented  at  the  1 6th  Explosive  Safety  Seminar, 
Hollywood  Beach,  Florida,  September  1974. 


PREFACE 

The  investigation  described  in  this  report  was  authorized  under  PA,  A 4932,  Project 
575 1 264.  The  work  was  performed  at  Southwest  Research  Institute  under  Contracts 
DAAD05-74-C-0751  and  DAAA15-75-C-0083. 

The  use  of  trade  names  in  this  report  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  of  such  commercial  hardware  or  software.  This  report  may  not  be  cited 
for  the  purposes  of  advertisements. 

The  information  in  this  document  has  been  cleared  for  release  to  the  general  public. 
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METHODS  OF  PREDICTING  BLAST  LOADS  INSIDE  AND 
BLAST  LOADS  OUTSIDE  SUPPRESSIVE  STRUCTURES 


I.  INTRODUCTION 

The  loading  froiii  an  explosive  charge  detonated  within  a vented  or  unvented  structure 
consists  of  two  almost  distinct  phases.  The  first  phase  is  that  of  reflected  blast  loading.  It 
consists  of  the  initial  high  pressure,  short  duration  reflected  wave,  plus  perhaps  several  later 
rctlectcd  pulses  arriving  at  times  closely  approximated  by  twice  the  average  time  of  arrival  at 
the  chamber  walls.  These  later  pulses  are  attenuated  in  amplitude  because  of  irreversible 
thermodynamic  process,  and  they  may  be  very  complex  in  waveform  because  of  the  com- 
plexity of  the  retlection  process  within  the  structure,  whether  vented  or  unvented.  If  the 
structure  has  solid  walls,  the  blast  loading  can  be  estimated  by  using  sources  of  compiled 
blast  data  lor  normally  retlected  blast  pressures  and  impulses  such  as  References  1 and  2.  and 
the  well-known  Hopkinson's  blast  scaling  law  (see  Chapter  3 of  Ref.  3).  The  effect  of  vented 
panels  in  the  suppressive  structures  on  reduction  of  the  reflected  blast  loading  can  be  very 
complex,  and  will  not  be  addressed  in  detail  in  this  paper. 

■•Vs  the  blast  waves  reflect  and  re-retlect  within  the  structure  and  as  unbumed  detonation 
products  combine  with  the  available  oxygen,*  a quasi-static  pressure  rise  occurs  and  the 
second  phase  of  loading  takes  place.  Proctor  and  Filler*  present  some  data  on  thse  pressures, 
Proctor*  has  developed  a computer  program  to  calculate  both  blast  and  quasi-static  pressure 
rises,  and  Sewell  and  Kinney*  also  present  methods  for  estimating  this  later  phase.  In 
addition,  Keenan  and  Tancreto’-  * have  made  measurements  of  blast  pressures  emitted  from 
rectangular  box  cubicles  with  various  vent  areas  and  pressure  rises  within  the  cubicles. 

Finally.  Lasseigne'*  has  measured  static  pressure  rises  in  closed  chambers  to  obtain  design 
information  for  a specific  suppressive  structure.  From  these  references,  one  obtains  the 
answer  that  for  the  particular  ratios  of  vent  area  to  chamber  volume  tested,  the  venting  has 
no  effect  on  the  peak  quasi-static  pressure.  Thus,  peak  static  pressures  for  unvented  or 
poorly  vented  structures  are  the  same.  Unfortunately,  essentially  no  data  exist  for  quasi- 
static pressures  within  well-vented  structures  and  the  crucial  question  of  the  actual  maxi- 
mum pressure  rise  within  such  chambers  remains  unanswered.  We  must  at  present  use  the 
unvented  pressure  rise  for  design  purposes.  We  have,  however,  conducted  a model  analysis 
and  fitted  curves  to  all  data  available  to  date  to  obtain  the  best  possible  estimate  of  this 
pressure.  The  model  analysis  and  curve  fits  are  presented  later  in  this  report. 

A third  important  question  regarding  blast  loading  and  suppressive  structures  is.  “Can 
blast  pressures  outside  these  structures  be  predicted  for  specific  designs?”  Many  of  the  past 
measurements  of  effectiveness  of  these  stnicturcs  have  been  based  on  blast  attenuation  which 
they  provide  (see  Refs.  9-1 2),  Using  these  references  and  more  recent  data  from  MTF.  we 
have  generated  a method  of  correlating  emitted  blast  waves  with  suppressive  structure  design 


'The  amouiu  of  oM^gen  avaiUbk  within  any  compleie  siruciuie  is  unatrccled  by  venting,  uniil  the  veniing  area  becumes 
very  large. 


7 


based  o.i  comparing  free-field  blast  data  to  blast  data  for  waves  emanating  from  supp.-essive 
structures.  This  method  introduces  an  effective  vent  area  ratio,  Oferf.  which  can  be  com- 
puted for  any  combination  of  vented  elements  in  a suppressive  structure  panel.  Using  this 
parameter  and  least-squares  curve  fits  to  free-field  and  suppressive  structures  blast  data,  we 
have  shown  that  the  influence  of  the  suppressive  structure  is  to  create  an  effective  standoff 
distance  R^,  . less  than  the  free-field  standoff  distance  Rj  at  which  side-cn  overpressure 
is  the  same  for  a given  blast  source  energy  W.  Alternatively,  this  method  will  predict  the 
i\‘duction  in  overpressure  over  a considerable  range  of  distances  outside  the  structure. 
Details  of  the  method  are  also  given  later. 


II,  BLAST  PRESSURES  OUTSIDE  SUPPRESSiVE  STRUCTURES 


The  side-on  overpressures  P^  in  the  free-field  around  an  explosive  charge  arc  given  by  a 
functional  relationship  as  expressed  in  Eq,  {I ). 


P.=f 


(free  field) 


(1) 


where 

R standoff  distance 

If'  - charge  weight 

Tliis  functional  relationship  is  the  famous  Hopkinson  blast  scaling  law  for  the  blast  field 
around  geometrically  similar  sources  at  sea-level  ambient  atmospheric  conditions,^  Assume 
that  a cubical  blast  suppressive  structure  whose  length  on  any  side  is  X and  whose  walls  are 
fabricated  of  a single  metal  sheet  with  holes  drilled  in  it  is  now  centered  over  the  explosive 
charge.  The  ratio  of  the  vent  area  of  a wall  to  the  total  presented  area  of  the  wall  will  be  de- 
fined as  equaling  a.  Equation  ( 1 ) for  free-field  blast  will  now  be  modified  by  the  additional 
geometric  parameters  defining  the  size  of  the  suppressive  cube  X and  the  vent  area  ratio  a. 

If  we  elect  to  write  a modified  form  for  Eq,  (1 ) in  nondimensional  ter.ns,  a functional 
equation  for  predicting  blast  pressures  outside  the  suppressive  structure  becomes: 


(suppressive  structure  equation) 


C) 


Equation  (2)  represents  a four-parameter  space  of  nondimcnsional  numbers  or  pi  terms. 
Although  no  functional  format  is  expressed  by  Eq.  (2),  sufficient  quantities  of  experimental 
data  can  be  used  to  obtain  an  empirical  relationship.  This  is  precisely  what  is  done  to  develop 
a relationship  for  predicting  blast  pressures  outside  of  the  suppressive  structure;  however, 
we  must  first  realize  that  most  suppressive  structures  do  not  have  walls  which  are  a single 
sheet  with  holes.  The  vast  majority  of  structures  have  three  to  six  wall  layers  with  various 
staggered  venting  patterns  so  fragments  will  not  escape  the  confinement.  This  means  that,  for 
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a multiwalled  confinement,  we  must  compute  an  effective  or.  ajff , so  Eq.  (2)  can  be  use  J to 
predict  blast  pressures.  To  compute  Oecr  for  a multiwalled  structure,  we  have  assumed  that; 


111  I 

— = — + — + + — 

a*fr  “i  “j  “.V 


wliere  V = number  of  elements  in  a suppressive  structure  panel.  Or, 


(3a) 


/•=/V 


- 

; = I 


(3b) 


Althougli  no  theoretical  proof  of  this  relationship  is  presently  possible,  it  does  reach  the 
appropriate  limits  for  small  and  large  numbers  of  plates.  For  example,  if  only  one  plate  is 
present.  Oeff  = a, . as  it  should.  If  an  infinite  number  of  plates  are  present,  o^ff  = 0,  with  the 
flow  completely  choked.  If  one  of  the  plates  is  solid  and  thus  has  a zero  a,  Oeff  = 0.  as  it 
should  If  all  plates  have  the  same  value  for  a,  Oeir  = oi/.V,  which  is  a number  smaller  than  a 
for  a single  plate,  as  would  be  expected.  In  each  member,  a is  defined  according  to  Eq.  (4) 


^ _ ^veni 


(4) 


For  plates,  the  meaning  of  this  definition  is  obvious;  however,  in  angles  and  louvres,  the  def- 
inition is  less  obvious.  Figure  1 defines  a in  a series  of  angles. 

In  a louvre,  we  use  a similar  definition  of  a,  except  that  the  o,  determined  on  the  basis  of 
Eq.  (4),  is  multiplied  by  a factor  equal  to  1/2.  This  factor  was  applied  because  the  data  of 
Reference  1 2 indicate  that  Icuvres  are  more  efficient  in  constricting  flow  than  are  plates  with 
holes.  Perhaps  this  is  explained  by  the  fact  that  the  entrance  of  a louvre  is  perpendicular  to 
the  entrance  of  a hole  in  a wall.  As  will  be  shown  later,  the  factor  of  1/2  appears  to  be  jus- 
tified by  a curve  fit  to  the  experimental  data.  Figure  2 illustrates  our  definition  of  a for  a 
louvre. 

Now  that  Ogff  has  been  defined,  we  are  prepared  to  develop  a functional  format  for 
Eq.  (2).  This  format  was  developed  by  assuming  that  Eq.  (2)  can  be  expressed  as 


If  logarithms  are  taken  of  both  sides  of  this  equation, 

(log  Pj } = (log  A)  + .V,  (log  + ;Vj  ^log  + .Vj  (log  o*ff) 


•Va 


(5) 


(()) 
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MGURI  : DEFINITION  OF  q IN  A LOUVRE 

The  equation  is  made  linear,  and  a least-squares  curve  fit  can  be  developed  by  stating  that; 


R X 

' '08  '°8  B- 


R 


1 flog  T 

•V. 
-V, 


^•J 


= Ilog/>J 


(7) 


Substituting  matrix  notation  yields: 

l^)l.Vl=|/’l  (8) 

and  a least  squares  curve  fit  results  for  log  A , .V, , A'j , and  A'j  or  the  ,V  matrix  when: 


[.-VI  = \L^\  [P\ 


(9) 


Experimental  test  data  from  References  9-12  were  used  to  develop  this  curve  fit.  The 
resulting  equation  is 


P,  = 97b, 3 


R^'^X'’^ 


(10) 
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wIktc 


/*j  side-oil  pressure  (psi) 

It'  eliarge  weight  (lb) 

R slaiidolT distance  from  charge  (ft) 
■V  width  of  suppressive  cube  (ft) 
a etfective  vented  area  ratio  (-) 


Figure  3 is  a plot  of  Eq.  ( 10)  versus  the  experimental  data  points  used  to  compose  this 
plot  Equation  ( 10)  appears  to  fit  the  test  results  excellently.  One  standard  deviation  for  the 
experimental  data  about  the  line  in  Figure  3 equals  18.6%,  which  is  only  slightly  worse  than 
would  be  obtained  for  free-ficid  data.  Because  this  is  a curve  fit  to  test  data,  Eq.  ( 10)  should 
only  be  used  when  input  conditions  fall  within  variations  in  individual  pi  terms.  The 
variations  included  in  test  results  were: 


0 0263  <Oe,f<  0.60 

0.323 < <1.77  (II) 

R 


4.27  ft/lb'  '^ 


< 17.5  ft/lb'  ' 


The  test  data  include  results  for  a wide  variety  of  panel  geometries  and  numbers  of  vented 
layers  in  each  panel.  These  range  from  as  few  as  two  layers  to  as  many  as  five,  and  com- 
binations of  spaced  angles,  zees,  perforated  plates,  and  louvres. 

It  is  interesting  to  curve  fit  free-Held  side-on  blast  pressure  data  from  References  ')-l  2 
using  the  same  procedure  over  the  same  range  as  for  the  suppressive  structure  blast  field  data 
The  resulting  equation  for  free-field  data  is 


= 976.3  (12, 

A comparison  between  Eq.  ( 1 2)  and  the  test  data  points  can  be  seen  in  F'igure  4.  The  stand- 
ard deviation  for  blast  pressures  in  the  free-field  is  13.1',?  which  is  only  slightly  better  than 
the  standard  deviation  for  the  suppressive  structure  blast  field  equation  Naturally,  Eq.  (12) 
should  only  be  applied  whenever  RjW''^  js  between  the  limits  established  by  Eq.  (II). 

If  one  compares  Eq.  (10)  for  suppressive  structures  to  Eq.  ( 12)  for  free-field  blast,  it  is 
immediately  apparent  that  the  influence  of  (he  suppressive  structure  is  to  create  an  effective 
standoff  distance  less  than  the  free-field  standoff  distance  at  which  blast  pressures  are  the 
same  for  a given  energy  release.  This  standoff  distance  with  a structure  suppressing  the 
blast  /?j,  is  related  to  the  free-field  standoff  distance  /?/  by: 


II 


limMftyHiriir'  '-■■'•T  ■'  ■ 


or 


Xin 


‘eff 


(13) 


(14) 


III.  PRESSURE  RISE  INSIDE  STRUCTURE 

In  this  section,  we  will  discuss  the  quasi-static  pressure  rise  within  a suppressive  structure. 
To  create  a solution,  we  will  first  perform  a model  analysis.  The  problem  is  envisioned  as  an 
instantaneous  energy  release  of  magnitude  W inside  a confined  volume  V.  A vent  area 
(Ocrr  A ) exists  through  which  internal  gases  can  escape.  We  are  interested  in  predicting  the 
internal  pressure  rise  p and  its  decay  as  functions  of  time  /.  Ambient  atmospheric  pressure 
Po  exists  initially  inside  and  outside  the  confined  volume.  To  define  an  equation  of  state  for 
the  gases  in  this  problem,  we  need  two  additional  parameters,  the  ratio  of  specific  heats  y 
and  speed  of  sound  <o-  Table  I summarizes  the  parameters  in  this  problem  and  lists  their 
fundamental  dimensions  in  an  engineering  system  of  force,  length  and  time  (F,  L,  T). 


Texts  such  as  Reference  1 3 tell  how  nond'meiisional  numbers  or  pi  terms  can  be  devel- 
oped from  this  list  of  variables.  Because  no  new  assumptions  are  inserted  in  developing  pi 
terms,  we  will  present  only  the  results  and  not  perform  all  of  the  algebra.  The  assumptions 
in  this  analysis  are  all  included  in  the  definition  of  the  problem,  so  that  phenomena  are  not 

considered  which  have  no  parameter 

TABLt  I PARAMETERS  DETERMINING  QUASt-STATIC 
PRESSURE  INSIDE  A VENTED  CONTAINMENT  VESSEL 


Ptramcier 

Symbol 

FuiKUmeniel 

Dimeiuioni 

Dcacription 

Volume 

V 

>Oaomctry 

Vented  Area 

<«efM) 

Energy  Rcleetc 

W 

FL 

Input  energy 

Atomoapheric 

PrcMure 

Po 

FIV 

1 

Sound  Speed  in 

1 

f 

Atr 

LIT 

^InlUel  tuie  of  tir 

Specific  Heii 

1 

1 

Rello  Alt 

7 

/ 

Pleat  ure 

IncrMie 

p 

FIL' 

• Retpotue 

Time 

1 

T 

> 

listed  in  Table  1 . Probably  the  major 
assumption  is  that  no  thermal  effects 
arc  considered;  in  other  words,  the 
pressures  dissipate  through  the  venting 
and  not  through  the  conduction  of  heat 
into  walls  of  the  suppressive  structure, 
therefore  an  acceptable  set  of  pi  terms 
which  can  result  is: 


JT.  = 


PaV 


JTs  =■ 


\n 


(15) 


If  we  assume  7 is  a constant  and 
are  only  interested  in  predicting  peak 
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pressure,  the  result  would  not  be  dependent  upon  time  or  the  pi  terms  ^3  and  jtj  would  be 
invariant.  Hence, 


Pmtx  ^ f / W 

■— — r, — ) 

Po  ^Po  y ' 

Because  \)q  is  also  invariant,  we  can  write  a dimensional  functional  format  for  Eq.  (16). 

Pmm  ~^i\y  ' y ) 

Figure  5 is  a plot  of  Phiix  versus  h'/F’  for  various  values  of  jtj  . Provided  Wj  is  less  than 
0.0775,  the  experimental  data  indicate  that  the  maximum  pressure  p^x  independent  of 
(tteff/l)’  VI  . This  can  be  written  as 

Pm.x  =/3  (p)  . < 0.077$  (18) 

The  data  used  in  developing  Figure  5 come  from  References  7 and  8.  In  addition  to  pre- 
senting their  own  data  which  were  obtained  at  the  Naval  Civil  Engineering  Laboratory, 
Keenan  and  Tancrcto  also  report  test  data  obtained  by  Proctor  at  the  Naval  Ordnance  Lab- 
oratory. Both  groups  of  experiments  used  Comp  B explosive,  but,  as  can  be  seen  in  Figure  5, 
their  experiments  were  in  different  domains  ot  W/V. 

The  dashed  straiglit  lines  in  Figure  5 are  the  asymptotes  for  complete  energy  conversion 
or  for /)„,,*  proportional  to  (WIV).  If  {WjV)  is  too  large,  insufficient  oxygen  is  available  to 
convert  all  the  energy  in  the  explosive  charge;  hence,  the  energy  release  is  reduced  by  the 
ratio  of  the  heat  of  detonation  divided  by  the  heat  of  combustion.  Figure  5 implies  that  for 
WIV  <,  0.003  complete  oxidation  occurs;  for  If / F > 0. 1 . the  only  oxidizer  available  is  that 
in  the  explosive  itself,  and  WIV  between  0.003  and  0.1  results  in  partial  afterburning. 

If  the  maximum  pressure  is  reached  before  significant  venting  occurs,  the  blow-down 
time  will  be  independent  of  ir4 , and  we  can  write  a functional  equation  for  time  of  blow- 
down. 


I2L  -f(P  \ 


(19) 


Neglecting  the  invariant  ambient  gas  parameters  in  Eq.  (19)  permits  us  to  write  a dimensional 
form  of  Eq.  (19). 
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FIGURE  5.  QUASI-STATIC  PRESSURE  RISE  INSIDE  AN  UNVENTED  ENCLOSURE 


The  data  used  to  develop  Figure  5 can  also  be  used  to  empirically  solve  Eq.  (20).  Figure  6 is 
a plot  of  ;/(/>'  '•  F'  ' ) versus  The  ordinate  of  this  graph  is  based  on  the  em- 
pirical observation  that  the  two  pi  terms  and  p can  be  combined  to  from 

We  can  now  write  Eq.  (20)  as  Eq.  (21): 


pUb  y\n  \ y } 

The  functional  format  for  Eq.  (21)  is  obtained  from  Figure  6. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  report,  we  have  presented  methods  for  estimating  long  term  (quasi-static) 
pressures  generated  by  internal  explosions  within  vented  suppressive  structures,  and  attenuated 
blast  pressures  escaping  from  suppressive  structures  with  various  vent  panel  design.  Scaled 
curves  for  appropriate  parameters  are  presented,  based  on  fits  to  available  experimental  data. 

We  have  several  suggest  ons  for  improving  methods  for  estimating  blast  loading  and 
attenuation  of  suppressive  structures: 

( 1 ) Initial  reflected  blast  loads  have  conservatively  been  estimated  to  be  the 
reflected  pressures  and  impulses  on  rigid,  non-vented  walls.  We  suggest 
that  shock-tube  or  field  tests  be  conducted  to  measure  these  loads  more 
accurately,  for  several  typical  vent  panel  designs. 

(2)  Experimental  data  for  quasi-static  pressure  rises  caused  by  internal  explosions 
in  vented  structures  have  been  limited  to  such  small  values  of  scaled  vent 
areas  that  the  maximum  pressure  rises  are  unaffected  by  the  vent  areas.  Tests 
should  be  run  on  well-vented  structures  to  determine  scaled  vent  areas  which 
which  cause  significant  reduction  in  quasi-static  pressures. 

(3)  As  additional  data  on  attenuated  blast  pressures  outside  suppressive  struc- 
tures are  obtained  in  the  course  of  subsequent  testing,  these  data  should  be 
factored  into  the  design  curves  and  equations  in  this  report  to  obtain  better 
fits.  In  addition,  curve  fits  should  be  made  to  scaled  data  for  blast  impulse 
outside  the  structures,  when  sufficient  data  are  available. 
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